Draft version February 2, 2008 

Preprint typeset using I^T^jX style cmulatcapj v. 6/22/04 



ANGULAR CORRELATIONS OF THE MEV COSMIC GAMMA RAY BACKGROUND 

Pengjie Zhang and John F. Beacom 

NASA/Formilab Astrophysics Center, Fermi National Accelerator Laboratory, Batavia, IL 60510-0500 

Draft version February 2, 2008 

ABSTRACT 

The measured cosmic gamma ray background (CGB) spectrum at MeV energies is in reasonable 
agreement with the predicted contribution from type la supernovae (SNIa). But the characteristic 
features in the SNIa gamma ray spectrum, weakened by integration over source redshifts, are hard 
to measure, and additionally the contributions from other sources in the MeV range are uncertain, 
so that the SNIa origin of the MeV CGB remains unproven. Since different CGB sources have 
different clustering properties and redshift distributions, by combining the CGB spectrum and angular 
correlation measurements, the contributions to the CGB could be identified and separated. The 
SNIa CGB large-scale structure follows that of galaxies. Its rms fluctuation at degree scales has a 
characteristic energy dependence, ranging from ~ 1% to order of unity and can be measured to several 
percent precision by proposed future satellites such as the Advanced Compton Telescope. With the 
identification of the SNIa contribution, the SNIa rate could be measured unambiguously as a function 
of redshift up to z ~ 1, by combining both the spectrum and angular correlation measurements, 
yielding new constraints on the star formation rate to even higher redshifts. Finally, we show that 
the gamma ray and neutrino backgrounds from supernovae should be closely connected, allowing an 
important consistency test from the measured data. Identification of the astrophysical contributions to 
the CGB would allow much greater sensitivity to an isotropic high-redshift CGB contribution arising 
in extra dimension or dark matter models. 

Subject headings: cosmology: large scale structure; gamma rays: theory-diffuse background; star: for- 
mation 



1. INTRODUCTION 

Type la supernovae (SNIa) produce intense fluxes 
of nuclear gamma rays, following the electron- 
capture decays of 56 Ni and 56 Co, contributing sig- 
nificantly to the MeV cosmic gamma ray back- 
ground (CGB). The predicted spectrum shape and 
normalization are both in reasonable agreement with 
CGB measurement s from COMPTEL l| Weidenspointnerj 

1999) and SMM l|Watanabe. Leising. Share, fc Kinzerl 

2000) . The case for the SNIa origin of the 
MeV CGB, including other possible sources, and 
the dependence on cosmological parameters and the 
star formation rate, has been nicely summ arized 
bv iWat anabe. Hartmann. Leising. fc Th el lll999| ) and 
iRuiz-LapuenteT'Casse'fcVa^^ . 

The SNIa gamma ray spectrum has strong line fea- 
tures around 1 MeV, and integrati on over SNIa redshifts 
turns these into a series of steps llClavton fc Silljll969l 
iClavton fc Ward! 119751 iThe. Leising. fc Clavtonl Il993j) . 
But insufficient signal to noise and residual systemat- 
ics of existing measurements smear these characteristic 
spectral features, and so the SNIa origin of the MeV 
CGB remains unproven. More fundamentally, combina- 
tions of several proposed contributions to the MeV CGB 
may also produce a spectrum in reasonable agreement 
with observations. Since redshift information is mixed 
when integrating over sources, this potential degeneracy 
makes identification of the CGB sources challenging. 

We propose new techniques to identify the origin of the 
MeV CGB, based on the observation that angular cor- 
relations of the CGB will reveal crucial information on 

Electronic address: zhangpj@fnal.gov, beacom@fnal.gov 



the three-dimensional source distributions. Since SNIa 
are associated with galaxies, the SNIa CGB should fol- 
low the galaxy distribution. Other possible M eV CGB 
sources, such as MeV blazars ijBJom et L _al.l ll995|) and cos- 
mic ray shock acceleration fMiniatil l20Q3? ) , have different 
clustering and redshift distributions, and would lead to 
different angular correlations in the CGB. MeV blazars, 
being rare and exhibiting strong variability ijZhang et alJ 
2002), will lead to uncorrelated angular fluctuations. The 
contribution from shock acceleration models would peak 
at low redshift (in contrast to the SNIa contribution 
which follows the star formation rate), and would not 
show features in the energy spectrum. Thus CGB an- 
gular correlations can be exploited to break degenera- 
cies that are present if only the flux is used. More im- 
portantly, the line emission features of the SNIa CGB 
are strongly enhanced in the autocorrelation function, 
allowing unambiguous identification of the SNIa contri- 
bution to the CGB. Furthermore, by cross correlating 
the CGB with galaxies, important redshift information 
can be recovered. The correlated angular fluctuations 
in the CGB that we consider arc intrinsically different 
from the uncorrelate d Poi sson fluctuations considered b y 
ICline fc Gaol ilOQCfl and IGao. Cline. fc Steckeil tlMl) : 
that case is appropriate if there are very few sources or 
detected photons per angular bin, which is not the case 
here. 

Such unambiguous extraction of SNIa CGB allows a 
statistically robust measure of the supernova rate (SNR) 
and its evolution and hence also the star formation 
rate (SFR). A routine one year CGB survey would re- 
flect the statistics of ~ 10 7 unresolved SNIa, with lit- 
tle obscuration by dust. In contrast, individual SNIa 
seen by gamma ray telescopes are expected to be very 
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Fig. 1. — SNIa spectrum, CGB mean flux spectrum F (in units of 
photons cm -2 s _1 sr — 1 keV — 1 ), angular power spectrum C 7 , and 
cross correlation power spectrum C 7 g with galaxies. C 7 and C 1 g 
are represented in two ways. Both scale as t~ 13 and are evaluated 
at I = 100. The CGB-galaxy cross correlation coefficient r is also 
shown. Emission line features are smoothed out in F(E) and C 1 g 
by integration over SNIa redshifts but survive in C 7 because the 
auto correlation amplifies these line features. The large amplitude 
of Cj at high E reflects strong CGB correlations at small scales 
and low redshifts. For the same reason, the CGB and galaxies have 
different redshift contributions, causing a decrement in C 1 g an d r. 



rare ijTimmes fe Wooslevl HqQ?). Even the proposed 
SNAP satellite would only optically resolve ~ 1000 SNIa 
per year. 

We use the SNIa contribution to the CGB to illustrate 
our proposed techniques. We calculate the SNIa CGB 
mean flux spectrum in the SNIa CGB auto correla- 
tion in ^3 and the cross correlation with galaxies in 21 
discussing the observational feasibility in S|5] Based on 
the measurements of flux spectrum, auto correlation, and 
cross correlation with galaxies, respectively, we present 
three new and independent methods to directly measure 
the SNR. Finally, we discuss the impact of improved sen- 
sitivity to new physics contributions to the CGB. 

2. CGB MEAN FLUX SPECTRUM 

The mean (angle-averaged) number flux spectrum of 
the SNIa CGB is 



F(E) 



L(E(l + z),z) 

47T 



-Rsn(z) 



fl 



where L(E) is the number of gamma rays per SNIa 
per energy interval, Rsn(z) is the SNR, and x is the 
comoving distance (we assume f2 m = 0.3, Q\ = 1 — f2 m , 
and h = 0.7) and dx/dz = {c / H Q ) / Q, m {l + z ) 3 + ^A- 
Since L(E) is wel l-predicted, F(E) can be used to con- 
strain the SFR dW atanabe, Hartmann, Leising, fc Thd 
119991 iRuiz-Lapuente. Casse fc Vangioni-Flaml [2001) . 
The SNIa spectrum, adopted from Fig. 4 of 



Fig. 2. — Logarithmic contributions (e.g., the integrands which 
yield the indicated totals when integrated with dz/z) to F (thin 
solid lines), C-yF 2 (thick solid lines) and C 7 g-F£<3 (dashed lines) 
from different z for the displayed energy values. These are domi- 
nated by line emissions in a series of redshift ranges, which allows 
inversion for the SNR. It also shows that C 7 g can be best mea- 
sured by choosing the redshift interval in the galaxy survey to best 
overlap the most significant SNIa redshift range. 



IRuiz-Lapuente. Casse fc Vangioni-Flaml l)200lD . is 
shown in the top panel of Fig. ^ The resulting mean 
flux spectrum, shown in the middle panel of Fig. ^ is 
sufficient to account for b oth the shape and no rmal- 
ization of the CO MPTEL llWeidenspoint ner|| 1999|) and 
(not shown) SMM (Watan abe. Leising. Share, fc KinzerJ 
2000) data over the energy range 0.5 < E < 3 MeV. At 
lower or higher energies, other sources besides SNIa are 
required to reproduce the observations. 

The SNR is determined by the SFR, the effi- 
ciency for producing SNIa, and the characteristic time 
scale tsN for a binary system to produce a SNIa 
(isN = 1 G yr for a white dwarf + non-white- 
dwar f system j|Ruiz-Lapuente. Casse fc Y angioni-Flarr 
1 200 11) 1. For the SFR. we adopt ijHippelein et al. 
12003(1 R S f{z < 1.2) oc exp(£/2.6 Gyr) and R S f{z > 
1.2) = R SF (1.2)exp(-(t-t(z = 1.2))/2.5 Gyr), both 
in units of Mq yr _1 Mpc~ 3 , and where t is 

the look-back time. Following the concl usions 

of iWatanabe. Hartmann. Leising. fc Thd (Il999f) and 
IRuiz-Lapuente. Casse fc Vangioni-Flaml l)2001tl that the 
SNIa contribution can likely account for the observed 
CGB, we assume this to be true and concentrate on 
normalization-independent techniques to test this hy- 
pothesis. We normalize the product of the star formation 
rate, efficiency to produce a SNIa, 56 Ni yield per super- 
nova, and gamma-ray escape probability to match the 
COMPTEL data. 

Past efforts of analyzing the mean flux spectrum can 
be improved by an inversion method. The spectrum can 
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be discretized as matrix product, 
F = K R SN 



(2) 



where the components of F and Rsn are F(Ei) and 
RsN(zj), and kernel matrix K is specified by 



K, 



dx/dz\ Zj 

47T 



Zj+Azj/2 



L(Ei(l + z),z)dz. (3) 



Zj —Azj/2 



Due to strong variations in L(E), we must use the inte- 
gral form instead of evaluating it at the median Zj . 

We split L into a smooth part L s and a line emission 
part L L . The line emission part L L can be approximated 
as a sum of delta-like functions, Yl L^(E, E^). The cen- 
tral values E^ do not change with z, but L^(E, E^ n ) may 
depend on z due to metallicity and environment evolu- 
tion. For simplicity, we neglect this possible redshift de- 
pendence, which can in principle be calibrated by SNIa 
observations in other bands. Thus 



K, 



dx/dz\ 

47T 



L s (Ei(l + Zj),Zj)Azj + 



E- 



Ah 



E, 



(4) 

where g m (z) — J L^(E , z) dE is the total number 
flux from the E^ n emission line at z, and the values m 
are determined by E^/E - 1 € (zj - Azj/2, zj + Azj/2). 
In the relevant range, F(E) is dominated by the line 
emissions, and the smooth part (L s ) can be neglected in 
Eq. 0. 

For each E, at most several specific emission lines and 
hence relatively narrow ranges of redshift contribute, as 
shown in Fig. [21 Since these peaks are narrow in redshift 
space, Azj — 0.1 is sufficient to minimize edge effects. 
It is sufficient to only keep the first several dominant 
g m , so that K reduces to a sparse matrix. This keeps 
the dominant signal while minimizing noise from other z 
ranges, simplifying and stabilizing the inversion. While 
present data are not sufficient, this method of inversion 
from F(E) can in principle be used to measure the SNR 
up to z ~ 1.5, as shown in Fig. O 

Data irregularities caused by statistical and systematic 
errors may make the inversion unstable, and more ad- 
vanced methods (e.g. iPen. Lu. van Waerbeke. fc Mellierl 
(2003)) or maximum likelihood analysis may help. We 
have assumed that the fraction of the CGB contributed 
by SNIa is near unity. If other sources contribute too, 
the SNIa fraction can be measured if the characteristic 
spectral features are seen. Near those features, the SNIa 
contribution will dominate, especially in the CGB auto 
angular correlation, as shown in Fig. Q and will be ex- 
plained in 

3. CGB ANGULAR AUTOCORRELATION 

The CGB angular autocorrelation has a different de- 
pendence on the SNR and provides a new method to mea- 
sure it. Since SNIa appear in all galaxy types, we expect 
the that the SNR follows the galaxy distribution with 
a constant bias b after averaging over a sufficient num- 
ber of galaxies. Additionally, since only SNIa produce 
the characteristic spectral features, the redshift averaged 
bias could be measured in principle. The angular reso- 
lution of the CGB experiments is generally ~ 1°, and at 
these scales we expect Gaussian fluctuations in the CGB. 



Since a (l ) 2 patch of sky contains ~ 10 5 galaxies, the 
constant bias assumption should hold very well. Then 
F(E) is given by 



F{E) = J 



L{E(1 



4tt 



-Rsn(z) 



(l+b8 G )d X , (5) 



where Sg is the relative galaxy number density fluctua- 
tion. For simplicity, we assume b — 1 throughout this 
paper. 

The widths of emission lines are generally several per- 
cent of the peak energies (see Fig. ^| , corresponding to 
spatial separations of several 30 Mpc/h. Since these are 
much larger than the galaxy- galaxy correlat ion length, 
we can use Limber's equation l)Peacockl ll999') in the rel- 
evant multipole range (t > 10 or < 30°) to obtain C 7 : 
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Due to the amplification in the square, only the line emis- 
sion features in L need to be considered, and g2,m is 
defined analogously to g m above, except for L 2 (E, z) in- 
stead of L(E, z). Since both the galaxy power spectrum 
(variance) A G and C 1 F 2 are observables, by the inver- 
sion methods discussed in |J2I Eq. © can be applied to 
obtain Rsn(z). Since L 2 is much sharper than L, the 
kernel matrix is more sparse than that of F, so that dif- 
ferent E ranges are less correlated, improving the error 
and the stability of the inversion. 

To quantify the sparse feature of this kernel matrix, we 
calculate the contributions to C 7 from various z. Here 
Aq can be calculated from galaxy correlation function, 
which we assume to be £,G( r , z ) — D 2 (z)(r/[5 Mpc/h]) a , 
with a = —1.7 , as inferred from the SDSS angular corre- 
lation results l|Connollv et aljr2002l) . and where D(z) is 
the linear density growth factor. These assumptions are 
accurate enough to illustrate our main results. 

A direct prediction is then C 7 (£) cx Cg(£) oc g~ a ~ 3 = 
l" 1 - 3 . As expected, C 7 shows strong signatures of the 
emission lines, as shown in Fig.^ There are only several 
important E^ (Fig. |5J), so the inversion is straightfor- 
ward. For example, a measurement at 0.5 MeV would 
directly measure Rsn(z) at z — 0.7. Since the density 
fluctuations weaken toward high z, C 7 is mainly sensitive 
to measuring Rsn(z) for z < 1. 

Here we briefly discuss the effects on our results due to 
variation of the parameters in our model. As noted, the 
angular power spectra are insensitive to the normaliza- 
tion of the CGB and hence the SFR, under the assump- 
tion that SNIa are the dominant contribution. Star for- 
mation histories weighted more heavily toward low red- 
shift produce larger angular power spectra, due to the 
growth of galaxy clustering. A larger delay time £sn 
causes SNIa to explode later, and would have a similar 
effect. Intrinsic variation in tsN would not affect extrac- 
tion of the SNR, but would introduce an additional un- 
certainty on the inferred SFR. The feasible redshift bins 
used in an inversion are unlikely to be narrower than 



4 



Az = 0.1, corresponding to a time interval of ~ 1 Gyr, 
so the error caused by fixing tsN should be negligible. 

The number of SNIa is finite, and so the measured 
correlations will reflect some Poisson noise, which can be 
subtracted. About N ~ 10 7 SNIa can be observed and 
thus the Poisson noise is about 1/yN < 10 ~ 3 , which 
is smaller than the angular correlation signal at degree 
scales (Fig. This estimation assumes equal weights 
in the Poisson noise for SNIa with different flux and at 
different z. In reality, nearby SNIa should be more heav- 
ily weighted, since their bright and rare occurrence will 
increase the Poisson noise. However, they can be sub- 
tracted as foreground sources. For a redshift bin cen- 
tered at z and spanning Az, in a one year survey with 
one steradian field of view, N ~ (l — 10) x 10 6 z 2 Az 
SNIa can be observed (the prefractor depends on the 
SNR at corresponding redshift), with fractional Poisson 
error ~ 10" 3 (z 2 Az)- 1 / 2 . The volume of a (1°) 2 angular 
bin is V ~ 8 x 10 6 (z 2 Az)(Mpc//i) 3 , and thus the density 
fluctuation in this volume is ~ (V 1 / 3 /8)~ 1 - 7 / 2 <t$D(z) ~ 
0.05(z 2 Az)~°- 3 D(z). Then, for the measurement of CGB 
flux from these angular and redshift bins, the S/N is 
~ 50(z 2 Az)°- 2 D(z). So, only for the most nearby red- 
shift bins with z < 0.01 and Az < 0.01, does the Poisson 
error exceed the signal. With a proposed gamma-ray 
telescope such as Advanced Compton Telescope, with 
sensitivity about 10~ 7 7 s _1 cm~ 2 , most bright nearby 
SNIa within ~ 100 Mpc (z < 0. 02) can be resolved 
and re moved, as can be inferred from lTimmes fc Wooslevl 
(1997). So, the Poisson error is at < 10% level. As we 
will see in this is not the dominant error source. So 
for the purpose of this paper, we omit the complexity in 
estimating this Poisson noise in more detail. 

4. CGB-GALAXY ANGULAR CROSS CORRELATION 

Since the SNIa CGB is tightly correlated with galaxies, 
by using the galaxy clustering data as a function of pho- 
tometric redshift range, the SNIa redshift distribution 
can be recovered. This cross correlation provides a new 
and robust way to identify CGB sources. For example, 
any CGB sources at high redshifts have no correlation 
with galaxies, which is particularly important for testing 
some new physics models, discussed further below. Fur- 
thermore, the cross correlation provides another way to 
measure the SNR. The cross correlation power spectrum 
CjG with galaxies in the redshift range [zi , Z2] is 
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where n(z) is the galaxy number distribution func- 
tion and the galaxy surface density is defined as 
Eg = / dn/dz{\ + Sc)dz. We adopt dn/dz = 
3z 2 /2/(z„,,/1.412) 3 exp(~(1.412z/z m ) 1 - 5 ) and choose the 
medi an redshift z m = 0.5 for SDSS (Do delson et al.1 
2002). 

Just as for F and C 7 , the contribution at a given energy 
is dominated by at most several redshift bins, as shown 
in Fig. |2 due to the SNIa line emission features. These 
redshift bins are then optimal for the cross correlation 
measurement because they contain the strongest signal. 



With large galaxy surveys, in which redshifts are mea- 
sured spectroscopically or photometrically (by colors), it 
would be easy to select the best ranges of redshifts for 
given energies. In the limit that the line emission features 
dominate, Eq. JJJ) becomes 

12 n <t\vv ^RsN(zj)dn gm{z 3 ) 
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which also allows a direct measurement of Rsn{z), where 
C 7 g is sensitive mainly to z < 1 because of the increased 
density fluctuations and more easily observable galaxies 
at lower redshifts. 

5. OBSERVATIONAL FEASIBILITY 

The three methods discussed above to recover FLgjy \z) 
involve very few assumptions. The cosmological param- 
eters have been or will be measured precisely, and mea- 
surements of galaxy clustering over a broad range of 
scales are also rapidly improving. The multi-wavelength 
properties and possible evolution of SNIa will be cali- 
brated by present and future observations. Thus the fea- 
sibility of our proposed methods depends mainly on the 
CGB observations. Low gamma ray fluxes and low S/N 
make the measurements of F(E), C^Qi and especially C 7 
highly challenging, but the results would provide new and 
important tests of i?sN {z) . Here we estimate the require- 
ments for such measurements. 

In the MeV range, the main noise sources arise 
from cosmic rays and rad ioactivities in the detec- 
tors l)Weidenspointnerj Il999|) ; these can be treated as 
Poisson noise and subtracted. Then the uncertainty in 
the measured C 7 is given by 

5{C 1 F 2 ) 
C 7 F 2 
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The noise power spectrum is given by C/v i7 — 
47r /sky B [( N N,-y/N 7 ) 2 /N N ^ + 1/Nj], where the first term 
is the Poisson noise of the instrumental background and 
the the second term is the Poisson noise of the signal. 
Here f2^ B l& tne fractional sky coverage, and iVjy i7 and 
A^ 7 are the total numbers of received instrumental events 
and CGB gammas, respectively. Here we have assumed 
equal weights for all SNIa. As discussed in fJH this sim- 
plification is sufficient for addressing N 1 , the most rele- 
vant parameter for the uncertainty on the CGB correla- 
tions. On the average, each SNIa contributes less than 
one detected gamma ray, so the signal Poisson noise is 
determined by A^ 7 instead of the total number of SNIa 
in the survey volume and observing period. The window 
function Wt reflects the angular resolution of CGB exper- 
iment, which we assume to be We = exp(— ^ 2 6 |2 /16 v / 7r), 
where 9 P is the angular resolution of the survey. 

For fixed survey time (fixed N^^ and AL), the optimal 
survey strategy has to compromise between a larger / s k y , 
which decreases the cosmic variance, and a smaller / s ky, 
which decreases the noise power spectrum Cat i7 . Thus 
there exists an optimal sky coverage f°£ for each £ such 
that the relative error of is the smallest: 
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The uncertainty in the measured C 7i g is then given by 
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The noise power spectrum of the galaxy survey is Cn, G = 
4nf^ y /N G , where N G is the number of galaxies observed 
with a sky coverage /^L, and r is the cross correlation 
coefficient between the CGB and a galaxy survey. For 
the cross correlation, larger f2^ B always improves the 
measurement. Our results are shown in Fig. |21 We have 
adopted SDSS as the target galaxy survey. We found 
that the optimal angular scale for correlation detections 
is the degree scale, essentially the smallest scale allowed 
by the detector angular resolution. Though higher ac- 
curacy could be achieved at smaller scales, some of the 
assumptions about Gaussian fluctuations and constant 
bias may begin to break down. At degree scales, for sur- 
veys containing N > 10 7 galaxies, the galaxy Poisson 
noise (the Cn, G /W^C g term in Eq. ITT1) is negligible, 

1/2 

and 5C lG /C lG cx / sky ■ Thus for a shallower survey 
covering the whole sky and containing the same amount 
of galaxies as SDSS, the accuracy of the measured C lG 
could be improved by a factor of 2. But in order to 
study high redshift SNIa with the cross correlation, sur- 
veys deeper than SDSS are required. 

The key requirements for the necessary detections are 
(a) good S/N (> 0.1), (b) large sky coverage (/ sky ~ 1), 
and (c) large number of gamma detections (iV 7 ~ 10 7 ). 
The gamma flux is ~ 10 7 (100 cm 2 ) -1 sr -1 yr _1 , so a 
field of view of ~ 1 sr and an integration time of ~ 1 
year are required. These requirements guarantee ~ 10 3 
contributing SNIa per pixel, and will average out the 
time dependence of the SNIa gamma ray emission (bright 
for ~ 1 year) and possible dispersion in SNIa spectra. 

COMPTEL satisfied all of these requirements but (a). 
Its poor S/N makes the measurement of C 7 and C lG un- 
feasible, as illustrated in Fig.0 The proposed Advanced 
Compto n Telescope (ACT) 1 satisfies all thre e require- 
ments {Milne. Kroeger. Kurfess. fc The! 120021) . With 
ACT and SDSS, several percent precision in both C 7 
and C 7 g measurement would be possible (Fig. [3J, and 
this would allow a c omparable precisio n in the inverted 
Rsn(z). MEGA 2 (IBloser et al.ll200l and and NCT 3 
fall in between COMPTEL and ACT. The ultimate pos- 
sible precision, ~ 1%, is set by the limited gamma flux 
in a one year survey. 

6. DISCUSSION AND SUMMARY 

We have presented the power of using angular corre- 
lations to identify and measure the SNIa component of 
the MeV CGB. The CGB fluctuations should follow from 
the underlying large-scale structure of galaxies on degree 
scales, with rms fluctuations varying from several ~ 1% 
to order unity, depending on the energy band. At the 
degree scale, such measurements are feasible and the in- 
terpretation is robust. We estimate that ACT + SDSS 
would be able to measure C 7 and C lG with several per- 
cent accuracy. The combination of CGB mean flux and 

1 ACT, |http://hese.nrl. navy.mil/gamma/det ector/ ACT/ A CT. htm| 

2 http:/ /www. mpe.mpg.de/gamma/instruments/mcga/www/mega 

3 http:/ /ssl. berkeley.edu/gamma/nct. html 
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Fig. 3. — Estimated optimal sky coverage (dashed line) and 
corresponding optimal relative errors in C 7 (solid lines) and C 1 g 
(dotted lines) measurements. We assume SDSS for the photo- 
metric galaxy survey, which will measure Nq ~ 5 x 10 7 galax- 
ies and cover a quarter of the sky. The SNIa CGB gamma flux is 
~ 10 7 (100 cm 2 ) -1 sr -1 yr" 1 . Since the power spectra vary slowly 
with I, we adopted a bin size A£ = 0.3£. Moving from left to right 
across the top row, and then again across the bottom row, the val- 
ues of (Nj,Nj/N NtJ ) are: (10 7 , 0.01), (10 7 , 0.1), (10 7 , 1), (10 s , 1), 
respectively in rough correspondence to COMPTEL, MEGA or 
NCT, ACT, and a several-year ACT-like experiment with better 
angular resolution. 



angular correlations will help to separate SNIa contribu- 
tions from other possible sources and allows an statis- 
tically robust measurement of the SNIa rate to z ~ 1. 
Additionally our proposed inversion technique using the 
mean flux spectrum would probe the SNIa and SFR rates 
to even higher redshifts. 

A variety of new physics models predict or al- 
low contributions to the present-day MeV CGB, 
for exam ple the dec ay of no n-baryonic cold dark 
matt er llOlive fc Sifkl 119851 iBarbieri fc Berezinskyl 
[1981 iDalvl 1198^ iGondolol Il99l lEllis et all [ 19921: 
" Kamionkowskil 119941: iKribs fc R.othsteinl 1199^ 
Abazajian^Fjfller^^F^atcl 2001: [Chen fc Kamionkowskil 



20041 iFeng. Raiaraman. fc Takavamall2003j) . the decay 
of massive gravit ons predicted by models of extra 

dimensions llArkani-Hamed. DimoDOiilos. fc Dvalil 

l999[ lHall fc Smithl 119991: lHannestad fc RaffeltJ 
2003; iCasse. Paul. Bertone. fc Siel 120041 



Kolb. Servant, fc I n ill 2003 1), and primordia l black 
hole evaporation llKim. L ee, fc MacGibb orJ 119991: 
[ Daghigh fc Kapustal 120021 iSendouda. Naeataki. fc Satol 
2003). When those photons arise from high-redshift 
sources, their angular distribution will be very isotropic 
and their energy spectrum may be featureless, though 
they may contribute an MeV CGB flux comparable 
to the observations. However, such sources may be 
htr&cparated from the SNIa contribution by using the tech- 
niques introduced in this paper, allowing more stringent 
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tests of new physics than with the flux constraint alone. 
Similarly, these techniques can discriminate against 
exotic gamma ray sources that follow the halo profile of 
the Milky Way. 

Finally, if the origin of the MeV cosmic gamma ray 
background is indeed type la supernovae, then this flux 
should be closely connected to the ~ 10 MeV cosmic 
neutrino background from type II, lb, and Ic (core- 
collapse) supernovae. Type la supernovae (arising from 
stars less massive than about 8M Q ) only efficiently pro- 
duce gamma rays, and core-collapse supernovae (arising 
from stars more massive than about 8M Q ) only efficiently 
produce neutrinos. Nevertheless, for an assumed initial 
mass function, the ratio of the gamma ray and neutrino 
fluxes should be mostly independent of the SFR. 

Each type la supernova produces ~ O.5M of 56 Ni, 
leading to ~ 10 55 gamma rays. In contrast, in a 
core-collapse supernova, the inner ~ 1.5-Mq core radi- 
ates about ~ 10 58 neutrinos, the dominant energy loss 
mechanism. The core-collapse supernova rate is about 
10 times larger than the type la supernova rate, so 
therefore the supernova neutrino background flux should 
exceed the supernova gamma ray background flux by 
~ 10 4 . Recently, the experimental limit on the su- 
pernova neutrino flux has been greatly improved by 
Malek et alJ (I2003f) . and is now very near the theoret- 
ical predictions jK aplinghat. Steigman. & Walken [2000 1: 
Fukugita & Kawasaki 2003:IAndo. Sato, fc Totani 2003; 



Strigari, KaDlingha^Sjejgman. fc Walkerll20f)l . In ad- 
dition. lBeacom fc Vaginsl ( |2004f) have shown that a mod- 
ification of the Super-Kamiokande detector would dra- 
matically improve the signal to noise, allowing a clear 
discovery. 

Though neither the gamma ray nor neutrino back- 
ground has been clearly measured yet, the predicted in- 
tegrated fluxes, with F v ~ 10 2 cm -2 s -1 (summed over 
flavors) and -F 7 ~ 10~ 2 cm~ 2 s _1 , are indeed in roughly 
the expected ratio. In the standard scenario, both back- 
grounds arise dominantly from supernovae with z < 1. 
Exotic scenarios might predict rather different flux ratios, 
or might produce gamma rays and neutrinos at high red- 
shift, such that the gamma rays cascade to lower energies 
by scattering, but the neutrinos simply redshift. Though 
at present this comparison is rather crude, it could be 
made much more precise, and the excellent prospects 
for detection of both the gamma ray and neutrino back- 
grounds makes this potentially a very important test of 
the assumed astrophysical scenarios (e.g., the timescale 
for a binary to form a SNIa) as well as more exotic physics 
models. 
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